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Non-gaussianity in the initial conditions of the universe is one of the most powerful mechanisms to 
discriminate among the competing theories of the early universe. Measurements using bispectrum 
of cosmic microwave background anisotropics are limited by the cosmic variance, i.e. available 
number of modes. Recent work has emphasized the possibility to probe non-gaussianity of local 
type using the scale dependence of large scale bias from highly biased tracers of large scale structure. 
However, this power spectrum method is also limited by cosmic variance, finite number of structures 
on the largest scales, and by the partial degeneracy with other cosmological parameters that can 
mimic the same effect. Here we propose an alternative method that solves both of these problems. 
It is based on the idea that on large scales halos are biased, but not stochastic, tracers of dark 
matter: by correlating a highly biased tracer of large scale structure against an unbiased tracer one 
eliminates the cosmic variance error, which can lead to a high signal to noise even from the structures 
comparable to the size of the survey. The error improvement on non-gaussianity parameter /„; 
relative to the power spectrum method scales as where P and n is the power spectrum 

and the number density of the biased tracer, respectively. For luminous red galaxies at z ~ 0.5 
we get an error improvement of about a factor of three over the current constraints, assuming 
unbiased tracers are also available with a high number density. For an ideal survey out to z ~ 2 
the corresponding error reduction can be as large as a factor of seven, which should guarantee a 
detection of non-gaussianity from an all sky survey of this type. The improvements could be even 
larger if two high density tracers with a significantly different sensitivity to non-gaussianity can be 
identified and measured over a large volume. 

PACS numbers: 98.80 



I. INTRODUCTION 



One of the major unsolved mysteries in cosmology is the creation of structure in the universe. There are many 
competing theories that differ in their predictions, some of which are accessible to cosmological observations today. 
Inflation pL, ^, 3, 4], the oldest and most successful of these theories, predicts that correlations of initial structures are 
nearly scale invariant, adiabatic, and nearly gaussian [H, 0, 0, [1, 0] • However, ever since the first detection of cosmic 
microwave background anisotropies in COBE it has been recognized that cosmic variance, a finite number of structures 
on the largest scales, provides a fundamental limitation to how well we can distinguish among competing models using 
two-point function correlations [Tol [Tlj . This is because the primordial density field is a single realization of a random 
process and only its two-point function is specified by the theoretical predictions. For many discriminating tests such 
as the shape of the power spectrum on large scales, the cosmic variance limit is fundamental. As a consequence, 
many interesting questions open today may remain unanswered, such as the question of whether the low quadrupole 
in CMB is due to manifestation of some new physics on very large scales or just a chance fluctuation 12]. 

Our desire to discriminate among the models motivates searching for other probes, specially those where the cosmic 
variance limit can be avoided. An important example where this limit can to some extent be avoided are primordial 
gravitational waves. This is due to the fact that there exists a polarization mode (so called B polarization) that only 
receives contributions from gravitational waves and not from the dominant scalar perturbations [3) [3 hj. In this 
case the limit is determined by instrumental and other irreducible noise rather than by the cosmic variance. Inflation 
allows for a detectable level of gravitational waves on the largest scales and as a result primordial gravitational waves 
have been declared as the smoking gun for inflation [l^ , based also on the fact that some competing models [l6l . ITtI . ITsj 
do not predict detectable levels of gravity wave signal and its detection would rule them out. The fact that cosmic 
variance can be avoided and that competing theories differ in their predictions are the primary motivators behind the 
interest in CMB polarization as the next experimental frontier. 

One would like to identify other probes that can discriminate among the competing theories. One of these is non- 
gaussianity and most of the models predict non-gaussianity of local type, ^ — ip + where $ is the gravitational 
potential in the matter era and (j) is the corresponding primordial gaussian case. In the simplest models of inflation 
primordial non-gaussianity is predicted to be very small, /„; 1. Alternatives to inflation, such as the new ekpy- 
rotic/cyclic scenario naturally predict large non-gaussianity, /„; ^ 1 [H, [13, ISI, leading to non-gaussianity as the 
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smoking gun for alternatives to the simplest models of inflation. However, all of the tests of non-gaussianity proposed 
so far suffer from cosmic variance limit. Until recently the most powerful method to place limits on /„; was based 
on bispectrum of cosmic microwave background (CMB), with the latest WMAP constraint —9 < /„; < 111 at 95% 
confidence level [2?| . The statistical power of this method is limited by the cosmic variance and with a better angular 
resolution one can sample more modes and thus improve the limits. For example, the limits should be improved 
significantly with the higher angular resolution Planck satellite, where 95% interval on /„; of 30 is expected [23], but 
this will be limited by astrophysical uncertainties such as the secondary anisotropics in CMB. Because the primordial 
CMB is damped on small scales it seems unlikely that the ultimate limits from CMB will be significantly improved 
over the current Planck predictions, but the promise of measuring high redshift 3-dimensional matter distribution 
with 21 cm transitions could improve the limits further, possibly to the level below unity ,24,]. 

An alternative approach using clustering of biased tracers of structure on very large scales has recently been 
proposed (25j . It was shown that non-gaussianity leads to a very unique scale dependence of the large scale bias, one 
that increases strongly towards large scales, and whose amplitude scales with the bias of the tracer relative to the 
dark matter. One can therefore place the limits on /„; by comparing the scale dependence of the power spectrum 
of the biased tracer to the one expected in cosmological models under the assumption of a scale independent bias. 
Subsequent work explored further theoretical issues and prospects for the future [2^ . A first application of 

this method, which we will call the power spectrum method, has recently been presented using large scale clustering 
of quasar and luminous red galaxies (LRG) galaxy data from Sloan Digital Sky Survey (SDSS) [13]. The result, 
—29 < /„; < 69 (95% C.I.), is already better than the latest CMB constraints from WMAP, suggesting this is a 
competitive method compared to the bispectrum from CMB and should be pursued further. Projections for the 
future suggest that a survey of 10^ LRGs out to 2: ~ 0.7 over a quarter of the sky, such as the proposed SDSS-III, will 
improve these limits by another factor of two [1^. Further into the future, an all sky sample of highly biased tracers 



out to z ^ 2 with 10^ galaxies, such as those contemplated for the dark energy mission [29], could reduce the width 
of the 95 % confidence interval to the level of 5-10 [28l[30l. [Slj. 

While these predictions for the future are impressive, it would be great if one could improve them further to reach 
the realm of a guaranteed detection, which is at the level of /„; ~ 1 [33] ■ There are two main obstacles to this. First, 
this method suffers significantly from cosmic variance: this is because the signal is strongest on the largest scales 
where the cosmic variance, due to the finite number of realizations at a given scale within a given volume, makes the 
error on the power spectrum large. Relative error from each mode is of order unity, hence we can only detect the fni 
signal if the relative change in power due to /„; is of order unity. Second limitation pointed out in j27|] is that for 
current data the effect of /„; is correlated with other cosmological parameters, such as matter density and primordial 
slope of the power spectrum, which also change the shape of the power spectrum on large scales. This is because 
over the limited range of scales accessible these parameters can partially mimic the effect of It is particularly 
important for two-dimensional projections, such as the quasar and LRG samples based on photometric data used in 
[27| . For a given multipole moment I the k-space projection window is broad and the effect of /„; extends to fairly 
high I, but it does not have a very characteristic scale dependence that would allow one to separate it from the other 
parameters. 

Above discussion suggests, firstly, that methods that can avoid the cosmic variance limit are particularly useful 
for future progress in discriminating between competing theories and, secondly, that non-gaussianity is a key probe 
that can differentiate among them. In this paper we propose a method to probe non-gaussianity from large scale 
structure that circumvents the cosmic variance limit and also eliminates the problem of its degeneracy with other 
cosmological parameters. Instead of measuring the power spectrum of a tracer and comparing it to predictions based 
on cosmological models we propose to compare directly the density field of a biased tracer to the one of an unbiased 
tracer (or, more generally, to a tracer with a different bias, which may include less biased or even anti-biased tracers). 
The relative bias of the two tracers is scale independent in the absence of non-gaussianity, but picks up a scale 
dependence in its presence. The main advantage is that tracers are generally biased, but not stochastic, on very large 
scales. This means that the precision with which the relative bias between the two tracers can be determined is only 
limited by the noise, which is given by the Poisson sampling of the field, and not by the cosmic variance. Moreover, 
because we are directly comparing two density fields any scale dependence of relative bias can only be caused by 
non-gaussianity and not by other cosmological parameters, which cannot affect the amplitude ratio even if they can 
affect the power spectrum. Hence, there is no degeneracy between /„; and other cosmological parameters with this 
method. In the next section we quantify these statements, computing the error predictions for both methods and 
comparing them to each other. We then apply it to some representative examples of present and future data sets, 
showing the improvements that can be expected. We finish by highlighting the requirements needed for the future 
surveys to maximize the power of this method. 
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II. ERROR ANALYSIS 



Let's assume we measure two tracers of matter density field, Si and 62- They arc both biased tracers of the underlying 
matter density field 6, Si = bid, where bi, 62 is the large scale bias of two tracers. We can introduce relative bias 
a = 61/62 and the corresponding covariance matrix elements in the Fourier domain are C22 =< S2 >= {P2 + ^2^) /V , 
C12 =< S1S2 >— raP2/V and Cn =< >= {a^P2 + nj~^)/V^, where r is the cross-correlation coefficient between 
the two fields, P2 is the power spectrum of second tracer and V is the volume element (which drops out from the final 
expressions) . 

We want to compare the errors on non-gaussianity parameter /„; as extracted from the power spectrum analysis to 
the one from the analysis of the relative amplitude of the two tracers. To answer this we need to compute the errors 
of the corresponding parameters, P2 and a, and their dependence on /„;. We will first compute the errors for these 
two parameters and then combine with their dependence on /„; to derive the final error predictions on /„;. Fisher 
matrix plays a key role in describing the ability of a survey to constrain parameters such as /„; . Its inverse gives the 
expected covariance matrix of parameters one wishes to estimate. It is defined as 

^^AA' = \tr[C.xC-^C,yC-\ (1) 

where C is the covariance matrix of the data defined above and A is the set of parameters one is estimating. 
Let us begin with the estimated error of the power spectrum P2 for a single mode. Applying equation [T] we find 

^ (0^X2 + Xif + 2(1 - r^)a^ [a\l - r^) + + + X2)] 
'''''' 2P|[a2(l-r2) + a2X2 + Xi+XiX2]2 ' 

where we defined Xi — {niP2)^^ ■ If Xi is larger than unity then we are in the Poisson noise dominated limit and 
Poisson noise dominates over the sampling variance. In the opposite limit <C 1 sampling variance dominates the 
error and the above expression reduces to Fp^p^ = = 2P| if r = 1, which is the usual sampling variance error 
for one mode. Hence the relative error from a single mode is limited from below to 2^/^ and this irreducible error is 
called the cosmic variance. 

We can also apply the above expressions to compute the diagonal component of the Fisher matrix for relative 
amplitude a, 

^ a^X2{l + 2X2) + r^Xijl + X2) + Q^(l - r^)(2 - + 3X2) 

[a^{l-r^)+a^X2 + Xi+XiX2]^ ' ^' 

which in the limit Xi <C 1, X2 ^1, 1 — ^ 1 becomes 

F-^ ^al^a^X2 + Xi+a^il-r^). (4) 

Thus the error on a from a single mode can be much less than unity if there is little stochasticity (r ^ 1) and the 
field is oversampled, Xi,X2 <C 1. Therefore, there is no fundamental cosmic variance limit here and if a depends on 
/„; then there is no cosmic variance limit on the latter either. 

To connect these expressions to the precision with which one can determine /„; we need to look at the dependence 
of P2 and a on /„; . For a given tracer with large scale bias 6 we have [2^ [2^ [l^, HI] 

P{k)^{b + Ab{k)fnl)^Pdm{k), (5) 

where Pdm{k) is the dark matter power spectrum and 

A6(fc) ^ 2{b-p)S^- = ^2k2Tik)D{z) ' 

where Sc = 1.686 is the spherical collapse linear over-density, (j) is the primordial potential in matter domination, S is 
the matter overdensity, Hq is the Hubble parameter, c is the speed of light, T(k) is the transfer function and D{z) is 
the linear growth rate normalized to (1 + z)^^ for z ^ 1. For tracers whose selection only depends on halo mass we 
have p — 1, but this can be larger if the tracer also depends on other properties such as the merging history, where 
for recent mergers one has p ~ 1.6 [2^ . Here we will use p = 1 in the following, but one should keep in mind that 
this is not always appropriate. 
Fisher matrix for /„; is 



(7) 
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Let us continue to treat the two ways of estimating /„; separately. From the power spectrum estimator for a single 
mode we find ^/„,/„,(^2) = Fp^p^{2Ah{k) /hP2)^ , which in the sampling variance limit gives 

Here we have assumed the non-gaussian correction is small to simplify the expressions. We have also assumed that the 
two samples have been combined optimally into a single sample with an overall bias h and non-gaussianity dependence 
A6, the details of which will depend on their number density and bias. For example, the simplest case is if 62 = 1 in 
which case we can assume this tracer contains no information on /„; and the overall sample is just the biased sample 
with b = bi and A6 = A61. The overall error is obtained by integrating over all modes 



''7'(p^) = ^l {^h{k)/hfedk, (9) 



where we assumed the sampling variance limit ignoring the Poisson term. Here V is the volume of the survey, 
kmin ^ 271 1 V'^^^ is the smallest wavector accessible in such a survey and kmax is the largest wavevector for which we 
can assume validity of this expression. Its value is often not very important since most of the sensitivity comes from 
large scales, ie small wavevectors. Note that we have ignored any degeneracies between and other cosmological 
parameters in this expression, so it is bound to be overly optimistic, although for anything more quantitative one 
would need to do the actual analysis with a given survey geometry and choosing the cosmological parameters one is 
varying in the analysis. 

The Fisher matrix for /„; from the relative amplitude of two tracers for a single mode is given by 

{a) = {X,^X,a-^ + l-rr'(^-^y . (10) 



Integrating over all the modes gives 



V 2 /^A&i(fc) A62(fc)^^ 



Here we have assumed that a = 61/62 is known perfectly from the smaller scales where /„; effects are negligible but 
linear bias still applies, so that information on a from large scales is used for the estimation of /„; rather than the 
relative bias itself. 

It is useful to compare the ratio of expected errors for a single mode in the sampling variance limit assuming r — \. 
We find 



A6/6 



A61/61 - A62/6 



(12) 



If we assume we have two tracers, one unbiased, 62 = 1, which contains no information on /„;, and one biased with 
hi = b this simplifies to 

^^4^ = ^2{X, + Xi/a^) = ^2[{n2P2)-' + (niPi)-i]. (13) 

The dependence on A6 and b has dropped out from this expression. Finally, if one assumes the abundance of 
unbiased tracers is much larger than that of biased tracers we can drop X2 term and we conclude that the signal to 
noise improvement per mode is simply given by the signal to noise power ratio of the biased tracer times root two, 
'^f,ii (^2) = (2/niPi)^^^, where Pi is the power spectrum of the biased tracer. 



III. DISCUSSION 



Let us apply the above derived relations to some examples of interest. SDSS-III plans to do a spectroscopic survey 
of about 1.5 X 10^ luminous red galaxies (LRGs) out to z ~ 0.7 over a quarter of the sky. The expected number 
density is about 3 x 10-''(h/Mpc)3 and average bias about b = 1.7 [111, so Pi{k = O.Olh/Mpc) = 0.8 x 10^(Mpc/h)3 
and we may expect about a factor of 3 error improvement if the unbiased tracers are available and their shot noise can 
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be neglected. This latter assumption is by no means easy to achieve, since lower bias galaxies live in lower mass halos, 
are typically fainter and SDSS-III will not target them spectroscopically. One way to pursue is using photometric 
samples, which is discussed further below. Another is to split the LRG sample into two luminosity bins, where we may 
expect for the brighter bin fii = 10-^(h/Mpc)3 and 5 ~ 2, so Pi{k = O.Olh/Mpc) 1.2 x 10^(Mpc/h)3 at z - 0.5, 
while for the fainter bin we may expect n2 = 2 x 10~'*(h/Mpc)'^ and b = 1.4, in which case da/dfni is reduced by more 
than a factor of two and the noise is increased since Xi = Xi = 0.1. In this case this method gives a comparable 
signal to noise to the power spectrum method. This may still be advantageous, since there is no degeneracy with 
other cosmological parameters and one can extend the analysis to small scales. 

Looking further into the future, there are a number of planned spectroscopic surveys that will measure a high 
number of redshifts to higher redshifts [13, H^- In this case one may expect larger improvements, specially if all of 
the most biased halos can be identified. For example, a,t z ^ 1.8 we may have ni ~ 3 x 10~'^(h/Mpc)'^ halos with 
average bias 6 = 2 and Pi{k = O.Olh/Mpc) ~ 3 x 10^(Mpc/h)'^ and most of these should be detectable with the above 
mentioned surveys. Hence in this case the improvement can be a factor of 7 over the power spectrum method. Note 
that with the power spectrum method such an improvement would require a 50-fold increase in the volume of the 
survey. This is a particularly exciting prospect since the predicted 95% confidence level interval for a full sky survey 
to z ~ 2 is of order A/„/ ~ 5 — 10 [28,, j50„ j51] and with the additional factor of 7 reduction of error we may be able 
achieve A/„/ ~ 1 (95% c.l. interval), at which point we enter into regime of a guaranteed detection [s^ . 

The benefits may be even greater for photometric surveys, where no spectroscopy is needed and hence can be done 
with a significantly smaller investment. On very large scales the advantages of a spectroscopic survey become less 
important, specially if the photometric redshift error is well below the largest scale of the survey and if, when several 
redshift slices are available, the cross-correlation information between photometric samples is included in analysis. 
Indeed, currently the strongest constraints come from the power spectrum method using LRG photometric sample and 
quasar photometric sample, both from SDSS j27| . As mentioned above, the main limitation of photometric samples 
is that the projection along the line of sight causes the non-gaussianity effect to extend over a wide range of multipole 
moments, making the effect difhcult to distinguish from other cosmological parameters using the power spectrum 
method [27]. This problem is eliminated using the method proposed here. 



Another potential advantage is that it is easier to obtain a high density sample of unbiased tracers at the same 
redshift using photometric data. For the LRG sample that goes to z ~ 0.7 the corresponding photometric low bias 
sample already exists from SDSS with well calibrated and reasonably accurate photometric redshifts and high number 
densities at the level of n = 10^ (number of galaxies per steradian) [33] . For photometric LRGs with their number 
density of n = 5 x 10^ and power C; = 1.5 x 10^'' at ^ = 6 [1^ we find a possible improvement of a factor of 5 over 
the current LRG photometric constraints of A/„i ^ 200 for 95% confidence interval. A potential challenge of this 
approach is matching the redshift distribution of the two samples, since any mismatch gives rise to less than perfect 
correlations. This can be achieved by constructing a set of weights for each galaxy that takes into account its redshift 
error distribution. It is also possible to combine spectroscopic and photometric samples. For example, one can use 
a high bias spectroscopic sample, such as the above mentioned LRG spectroscopic survey from SDSS-III, and apply 
the radial weights to reproduce the redshift distribution of a photometric sample. Finally, we should also note that 
the ultimate unbiased tracer insensitive to /„; to correlate against can be the dark matter itself, as measured from 
weak lensing, although in this case the radial window is rather broad [ssj . Future photometric surveys may allow for 
much more significant improvements. For example, PAU survey |39( plans to measure galaxies with a number density 
of n = 10^"^(h/Mpc)'^ out to z ^ 0.9, so here both the biased and unbiased samples will be available. The expected 
accuracy of photometric redshifts, dz ^ 0.003(1 + z), is sufficiently small that it can be ignored for the large scales of 
interest here. The expected error improvement will be about a factor of 3 over the power spectrum method and the 
predicted 95% c.l. interval of Afni ~ 5 is approaching the levels of a guaranteed detection. 

These predictions are based on the assumption that there is no stochasticity in the tracers, ie we assumed r = 1 in 
the analysis. It is well known that halos and galaxies are not perfect tracers of the density field and on smaller scales 
stochasticity develops. On small scales and/or for high density samples this can become the limiting factor: as shown 
in equation 1131 this happens when 1 — exceeds {nP)~^ and since we need the latter to be well below unity for the 
method to be superior over the power spectrum method this is a potential issue that needs further investigation. The 
large scale cross-correlation coefficient between different tracers has not been studied in much detail with numerical 
simulations. In [4^] the related cross-correlation between final dark matter and halos has been investigated, but 
the small box size used in the analysis makes the results noisy. Still, the results suggest that the cross-correlation 
coefficient is indeed small on very large scales, with 1 — < 0.1 for k < 0.05h/Mpc and likely to be even smaller 
on larger scales. An analysis with a much larger volume is in progress and preliminary results in configuration space 
suggest that 1 - ~ 0.1 at 5 Mpc/h and 1 - - 0.01 for r > 50 Mpc/h [4l|. This issue clearly requires further 
attention and more detailed studies using very large volume simulations with biased tracers will be needed to address 
it in detail to see if cross-correlation coefficient is sufficiently close to unity on very large scales to allow the several 
fold reduction of error promised by this method. 
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While we have focused on the biased tracers of large scale structure, it is also possible to look for the signature 
using the anti-biased tracers. For low mass halos the bias is asymptotically approaching b ~ 0.7 [13, S^l and for this 
population the effect of /„; has the opposite sign from the biased population with an fni dependence that is about 
a factor of 3-4 smaller than for b — 2 population assuming p — I. In addition, such an anti-biased tracer has about 
a factor of 10 smaller amount of power. These two disadvantages can in principle be offset by the higher number 
density, which can be orders of magnitude larger for b ^ 0.7 population relative to the b ^ 2 population. We can also 
use these antibiased tracers instead of the unbiased tracers and comparing them to the biased tracers, which improves 
the sensitivity to /„; and reduces the shot noise relative to the unbiased tracer we assumed in our analysis. 

More generally, recent studies suggest that bias can depend on quantities other than the halo mass [H, [43| , and 
similarly one can also expect secondary parameters that may enhance or suppress the sensitivity to /„;. Indeed, an 
extended Press-Schechter analysis suggests one such second parameter that suppresses the sensitivity to /„; of highly 
biased tracers is the recent merger activity [13] . Observationally identifiable candidates that satisfy this criterion may 
for example be quasars. The challenge for the future is to find two tracers, one that exhibits a significant dependence 
on fni and one that does not, both of which come with a sufficiently high number density to reduce the Poisson noise 
and both of which can be observationally identified within existing and future surveys. Given the potentially huge 
payoff of finding such tracers for the questions of interest to the fundamental theories of the universe this is a challenge 
that is worth exploring further both with observations as well as with numerical simulations and analytic methods. 

I thank Nikhil Padmanabhan and Anzc Slosar for useful comments on the draft. U.S. is supported by the Packard 
Foundation and Swiss National Foundation under contract 200021-116696/1. 
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